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Recently a class of supersymmetric gauge theories have been successfully implemented on the 
lattice. However, there has been an ongoing debate on whether lattice versions of some of these 
theories suffer from a sign problem, with independent simulations for the ,yV — (2,2) supersym- 
metric Yang-Mills theories in two dimensions yielding seemingly contradictory results. Here, 
we address this issue from an interesting theoretical point of view. We conjecture that the sign 
problem observed in some of the simulations is related to the so called Neuberger 0/0 problem, 
which arises in ordinary non-supersymmetric lattice gauge theories, and prevents the realization 
of Becchi-Rouet-Stora-Tyutin symmetry on the lattice. After discussing why we expect a sign 
problem in certain classes of supersymmetric lattice gauge theories far from the continuum limit, 
we argue that these problems can be evaded by use of a non-compact parametrization of the gauge 
link fields. 
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Introduction: There exist at least two major reasons to study supersymmetric gauge theories 
regularized on a space-time lattice. Firstly, supersymmetric Yang-Mills (SYM) theories, an im- 
portant subclass of these gauge theories, play prominant role in the AdS/CFT correspondence in 
which it is conjectured that string theories in AdS backgrounds can be described via a conformal 
supersymmetric gauge theory living on the boundary of the AdS space. At finite temperature this 
holographic duality relates properties of black holes in the AdS space to the thermal behavior of 
this holographic dual gauge theory. This duality is such that the supergravity limit of the string 
theory corresponds to the strongly coupled regime of the gauge theory. In such a situation one can 
hope to study the non-perturbative regime of the gauge theory using lattice methods - see for 
some recent lattice results. Secondly, supersymmetric field theories are among the most promising 
candidate models for the physics beyond the Standard Model but many questions related to su- 
persymmetry breaking are inherently non-perturbative in character and call for study using strong 
coupling methods such as lattice simulations. 

In spite of these obvious motivations, httle progress has been made in studying supersymmetric 
theories on the lattice until relatively recently. The problem is that naive discretizations of super- 
symmetric theories break supersymmetry completely leading to a fine tuning problem to regain 
supersymmetry as the continuum limit is taken. However, recently a number of new approaches 
have been advocated that address this problem and new classes of supersymmetric lattice gauge 
theories constructed - see [^, |3|, ^ ^, ^, ^ and [|8[ ^, 10|. In this paper we focus on just one of these 



new lattice formulations in which a part of the algebra of (extended) supersymmetry is maintained 
under discretization. This approach can be used to directly formulate a lattice theory with = 4 
supersymmetry in four dimensions and can be derived in two seemingly different ways. The first 
proceeds via orbifold projection of a matrix model iQ] while the second approach employs twisting 
of the spinorial supersymmetries of the target (continuum) SYM theory [||, ^ ^] . However it has 
been shown that both yield equivalent lattice constructions [11, 12, 13, |14[ ^. The lattice theories 



constructed this way are gauge invariant, preserve a subset of the original supersymmetries and 
provide the target (continuum) theories in the limit of vanishing lattice spacing. These construc- 
tions are successful only in implementing certain SYM theories with extended supersymmetries 
that include an extensive list of interesting class of SYM theories such as c/K = 2 in 2D, c/K = 4 in 
3D and 4D, and their dimensional reductions. 

However, even after implementing the SYM models on the lattice we potentially encounter an 
additional, and rather crucial, difficulty - the fermionic sign problem: Consider the lattice theory 
with a set of bosons and fermions (1/7, i/a). The partition function of the theory is 

Z = j [d4][dW][dw]^xp(^-SBm -WM[<^]\lf) = j mPf{M) exp (-5b[0]) , (1) 

where M is the fermion matrix which is an anti-symmetric matrix and Pf{M) is the Pfaffian of 
M. For a 2?i X 2?i matrix M, the Pfaffian is explicitly given as Pf{Mf = DetM. The Pfaffian, 
being a complex quantity, has a phase angle a that may change for different configurations. The 
a-weighted average then may cancel between different configurations resulting in a huge statistical 
error. This is the famous (fermionic) sign problem. 

In the lattice supersymmetry constructions, there has been an ongoing debate on the existence 
of the sign problem in the two-dimensional four supercharge {£2 = 4) lattice SYM theory. The 
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resolution of this sign problem is very crucial as the extraction of continuum physics from this 
lattice model depends very much on it. In [15], it was shown that there did exist a sign problem in 
the lattice J2 = 4, D = 2 SYM. However, since those results were obtained without any importance 
sampling of the partition function they would not necessarily hold in the full theory. However, 
in JTB]] it was shown that for the phase quenched ensemble at non zero lattice spacing there still 



appeared to be a sign problem for this theory. Later, Hanada et al. [17| showed that there was no 
sign problem for this theory in the continuum limit. In this work, we argue that for this theory, one 
must expect a sign problem in this model due to an extension of the so called Neuberger 0/0 prob- 
lem if one employs the usual (group based) parametrization of the gauge fields. However, if one 
employs a different non-compact parametrization of the bosonic fields, such as those used in [17] it 
may be avoided. In the following, we first explain why the usual gauge-fixing Faddeev-Popov (FP) 
procedure fails in lattice field theory giving rise to the conventional Neuberger 0/0 problem. Then 
we show that the partition function for the J2 = 4 SYM theories, being ^-exact, can be viewed as 
being obtained through a Faddeev-Popov gauge-fixing procedure and hence potentially subject to 
the Neuberger 0/0 problem, the latter appearing as a fermionic sign problem. However, we argue 
that the usual Neuberger problem only occurs if the relevant group manifold contains non-trivial 
cycles - a feature which is only true if an exponential parametrization of the gauge links is made; a 
non-compact parameterization side steps the issue. 

The Faddeev-Popov procedure and the Neuberger 0/0 problem on the lattice: Let us 
consider a gauge theory on a fiat Euclidean space-time lattice defined by the partition function 
Z = J[DU]exp{—S[U]), where the gauge fields Ua{n), defined along the links (n,n + /i^), with Ji^ 
being the basis vectors on the lattice, take values from a group G. Since the action S[U] possesses 
a gauge symmetry Ua{n) — > g{n)Ua{n)g^ (n + lla), there are infinitely many degrees of freedom 
for the gauge configurations and hence the partition function is ill-defined. Here, all the gauge 
transforms g{n) also take values from G. To remove the redundant degrees of freedom, one has 
to choose a representative out of the infinitely many possible physically equivalent configurations 
for a given set of gauge fields Ua{n)- This can be achieved by imposing a constraint on the gauge 
fields, called a gauge-fixing condition. In the perturbative limit, there is a well-defined procedure, 
the FP procedure, to impose the gauge-fixing condition in the path integral: to impose a condition 
f{U,g) = 0, one can insert a unity defined by Zgf •= 1 = / Dg5{f{U,g))Afp, and then integrate 
out the gauge transforms. Here, App is the Jacobian determinant of f{U,g) with respect to the g 
variables, and Zgf is called the gauge-fixing device. Writing Afp in terms of Grassmann variables 
(ghost and antighost fields), we obtain a gauge-fixed action which possesses a symmetry among 
gauge fields and ghosts and anti-ghosts, called the Becchi-Rouet-Stora-Tyutin (BRST) symmetry. 
Letting ^ be the BRST operator, the gauge-fixed part of the action can be written as Sgf = ^A, 
where A is a combination of terms involving gauge fields, ghosts and anti-ghosts. This feature of 
Sgf is called .^-exactness (and Sgf is called a J2-exact term). 

A crucial assumption here is that f{U,g) = should have a unique solution. While in the per- 
turbative limit, because f{U,g) = turns out to be a linear equation, there is a unique solution for 
such a condition, Gribov showed that in the non-perturbative regime for any covariant gauge-fixing 
condition there is in general more than one solution, the additional solutions are called Gribov 
copies 1 18]. 

On the lattice with standard, i.e., group based Wilson parametrization of the gauge fields. 
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there is another problem: the Gribov copies come in pairs with opposite signs of A^p. Thus the 
gauge-fixing partition function when summed over all Gribov copies is zero: Zgf = 0! Then, the 
expectation value of any gauge-fixed observable turns out to be 0/0; this is known as the Neuberger 
0/0 problem [19]. This problem has a very important consequence: a BRST symmetry can not 
be realized on the lattice for the standard Wilson formulation. This problem was interpreted by 
Schaden in terms of a Witten-type topological field theory [ ^0| ] and he arrived at the conclusion 
that on the lattice Zgf computes a topological invariant, the Euler characteristic, of the group 
manifold at each lattice-site, which is zero for SU (N) or U (N) (Also see [22, 23] for a recent 
point of view and [ 24 ] for the continuum counterpart.) 

The lattice SYM partition function as a gauge-fixing partition function: The 2D J2 = 4 
SYM theory is also a cS-exact theory, i.e., the complete action of this theory can be written as J2 
acting on some combination of bosons and fermions. In fact, this is a characteristic of a Witten- 
type topological field theory [2C]. In this case the operator J2 (the nilpotent BRST charge of the 
topological field theory) arises via twisting of the continuum supersymmetries. In other words, 
we can construct the complete two-dimensional J3 = 4 SYM theory described in the twist and 
orbifold prescriptions via gauge-fixing the fields to fiat complexified connections. We 

start with a trivial classical action Sdafisicaii-s^a, a] = 0, which exhibits a topological shift sym- 
metry S!/a — > -B^a = -^0 + Oa with da being an infinitesimal (complexified) shift. The topological 
symmetry is then gauge-fixed by following the Batalin-Vilkovisky [ p^ ] procedure. This process 
will generate a set of ghost, ghost-for-ghost and anti-ghost fields as end products, providing the 
necessary bosonic and fermionic terms in the action. We then relax the restrictions on the nilpotent 
BRST charge J2 and the topological nature of correlation functions to enhance the topological field 
theory to a supersymmetric gauge theory. The final action we achieve through this process is the 
supersymmetric lattice action. Thus the corresponding partition function of the twisted supesym- 
metric lattice theory can be viewed as the corresponding Zgf in the FP procedure for the classical 
action ^classical = 0. Another important feature of the Witten-type topological field theories is the 
existence of a Nicolai map: after integrating out the fermionic fields in the partition function (i.e., 
the one given in Eq. (|I|)), for the Witten-type topological field theories, there exist a change of 
variables (i.e., a map) that trivializes the bosonic action Sb[<P] and the Jacobian of that change of 
variables cancels the Pfaffian up to its sign. The final expression of Z can be shown to compute 
a topological invariant, namely the Euler characteristic (Since in the twisted SUSY theories, the 
gauge-fields are complexified, the actual topological invariant the partition function computes may 
not be precisely the Euler characteristic but some closely related topological quantity.). This means 
that this theory, when put on the lattice, should suffer from (an extended version of) the Neuberger 
0/0 problem. The difference here is that the Gribov copies now correspond to vacuum states of the 
supersymmetric lattice theory. The Neuberger 0/0 problem rephrased a la Schaden would imply 
that these vacuum states come in pairs with opposite signs of the Pfaffian and cancel each other - a 
sign problem. 

Differences between the two contemporary simulations approaches: The preceding argu- 
ment assumes that the lattice gauge fields are represented as exponentials of the continuum fields. 
Rather remarkably, these supersymmetric lattice constructions allow for another possibility; a non- 
compact parameterization based on expanding the fields on the algebra of the corresponding U (N) 
group. In the first approach the complexified gauge fields appearing in the twisted continuum the- 
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ory s^a{^) = ^a{x) + iBa{x), are mapped to link fields '^a(n) living on the link (ii,n + /i^,) through 
the mapping %,{n) =e<(n). The gauge links belong to the group G, which can be taken as the 
complexified U (N) group. We call this parametrization the exponential parametrization. However, 
since the fields live on complexified U (N) another possibility is open to us; namely we can simply 
expand the fields on the algebra of U {N) and require that the theory generate a vacuum expectation 
value for the imaginary part of the trace mode of the link fields. That is, '^^(n) = 1 + jS^(n). In 
order to not break U (N) gauge invariance one must then supplement the action with a gauge in- 
variant potential which automatically generates the required vev for the trace mode (such a vacuum 
expectation value is compatible with the moduli space of the theory). Since this expansion agrees 
with the usual expansion of the exponential in the Wilson formulation the two parametrizations 
lead to the same naive continuum limit. However, they can clearly be very different at non-zero 
lattice spacing. The exponential representation was used in [[I^, while in the non-compact 
representation was used. 

We have added a potential term of the form Sm = 1^'^ i:„[(^)Tr('^/(n)^a(n)) - 1]^ to the 
lattice action, with /i a tunable mass parameter, which can be used to control the expectation values 
and fluctuations of the lattice fields. Notice that such a potential obviously breaks supersymmetry 
- however because of the exact supersymmetry at = all supersymmetry breaking counterterms 
induced via quantum effects will possess couplings that vanish as /X — )• and so can be removed by 
sending /i — )• at the end of the calculation. 

We conjecture that the sign problem for the two-dimensional twisted SYM with exponen- 
tial parametrization is nothing but an extended Neuberger 0/0 problem. For the non-compact 
parametrization we live only in the algebra and hence the topology of the group manifold is irrel- 
evant. Interestingly, independent of the lattice SYM scenario, non-compact gauge fields have also 
recently been used to restore BRST symmetry on the lattice [22]. There the decompactification of 
the gauge fields is done via stereographic projection of the original compact gauge-fields. It should 
be noted that the Neuberger 0/0 problem applies strictly only to the lattice theory and it is not clear 
how it behaves as we approach the continuum limit. Note that the Neuberger 0/0 problem is a 
sufficient condition to have a sign problem in the corresponding theory and not a necessary condi- 
tion, i.e., there may still be sign problem in theories where the extended Neuberger 0/0 problem 
does not apply. In the case of two- and four-dimensional J2 = 16 twisted SYM, the action has 
a =S-closed piece in addition to a =S-exact piece so Neuberger will not apply independent of the 
parametrization of the gauge links. Of course, this theory may still have a sign problem originating 
from some other source although recent results in fact indicate that this theory is safe from a sign 
problem for both zero and non-zero lattice spacing |26]. Below we provide some results from the 
recent numerical simulations to support our conjecture. 

Simulation results: To simulate the theory, we rescale all the lattice fields by powers of 
the lattice spacing to make them dimensionless. This leads to an overall dimensionless coupling 
parameter of the form N/ {2Xa^), where a = p/T is the lattice spacing, j8 is the physical extent of 
the lattice in the Euclidean time direction and T is the number of lattice sites in the time direction. 
The lattice coupling is defined as fc = {NT'^) / {IX^'^), for the symmetric two-dimensional lattice, 
i.e., the spatial length L = T . Note that A is the dimensionless physical 't Hooft coupling. In our 
simulations, the continuum limit can be approached by fixing Aj8^ with N fixed and increasing the 
lattice size L — > oo. We take our parameters to be A = 1 .0, 2.0 and L = 2, • • • ,12. 
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Figure 1: The histograms ctf the phase angles a for ^ — A theories with gauge group t/(2),A = 1.0,;U=0.1, 
for non-compact representation. 
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Figure 2: The histograms of the phase angles a for = 4 theories with gauge group U{2), X — 2.0, /i = 0. 1, 
for exponential representation. 

We examined the Pfaffian phase measured in the phase quenched ensemble generated over 
sequential runs of about 20000 Monte Carlo configuration (MCC) time steps for each L and A. 
Measurements were taken after every 5 MC time steps. The simulations we show were performed 
for anti-periodic (thermal) boundary conditions for the fermions, and for both non-compact and 
exponential parameterizations of the gauge links. The first plot shows the expectation value of the 
absolute value of sin a (a is the Pfaffian phase) for the J2 = 4 model at = 0.1 using non-compact 
links together with a histogram of the phase angle. Clearly no sign problem is visible even on very 
coarse lattices. The second plot shows the same quantities using exponential links and the sign 
problem is visibly severe (in fact at its worst) away from the continuum in this case. 

Conclusions: In this work, we have discussed the possibility of a sign problem in two- 
dimensional lattice SUSY with twisted supersymmetries. It appears that the severity of any sign 
problem is heavily dependent on the parametrization of the gauge fields. For the non-compact 
parametrization, we show that there is no sign problem even on the smallest lattices while for the 
exponential parameterization the theory suffers from a severe sign problem away from the contin- 
uum limit. We have provided a possible theoretical explanation of this observed sign problem: with 
exponential parameterization, it is a manifestation of the Neuberger 0/0 problem. 
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